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Abstract: The spectral characteristics of a fiber Bragg grating (FBG) 
with a transversely inhomogeneous refractive index profile, differs con- 
siderably from that of a transversely uniform one. Transmission spectra 
of inhomogeneous and asymmetric FBGs that have been inscribed with 
focused ultrashort pulses with the so-called point-by-point technique are 
investigated. The cladding mode resonances of such FBGs can span a full 
octave in the spectrum and are very pronounced (deeper than 20dB). Using 
a coupled-mode approach, we compute the strength of resonant coupling 
and find that coupling into cladding modes of higher azimuthal order is very 
sensitive to the position of the modification in the core. Exploiting these 
properties allows precise control of such reflections and may lead to many 
new sensing applications. 
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1. Introduction 

In a typical transmission spectrum of fiber Bragg gratings (FBGs), it is common to observe a 
comb of resonances on the shorter wavelength side of the Bragg peak. These dips result from 
the resonant coupling between the forward-propagating core mode and modes that propagate 
backwards within the cladding of the fiber [1,2]. Since modes guided in the cladding of a fiber 
are very sensitive to bending of the fiber as well as the surrounding environment, they are widely 
exploited for sensing applications. In such sensor arrangements, FBGs and long period gratings 
(LPGs) are used as core to cladding mode converters [3]. However, conventional FBGs, in 
which the cross-section of the core is uniformly modified, couple light only into cladding modes 
of the lowest azimuthal order [1]. In general, coupling into higher azimuthal fiber modes can 
be realized by FBGs that are tilted [4-6] or have an asymmetric transverse index profile [7]. 
More elaborate in-fiber devices to achieve strong cladding mode resonances involve tapered 
fiber sections [8] or deliberately misaligned fiber splices [9]. 

In this paper, we investigate how strong cladding mode coupling can be achieved and con- 
trolled with highly localized FBGs. Such FBGs can be fabricated using focused ultrashort lasers 
operating in the near infrared (NIR) [10]. Due to the nonlinear nature of the absorption process 
of light pulses in the femtosecond regime, highly localized refractive index profiles can be ob- 
tained in non-photosensitive materials without damaging the surrounding fiber material. The 
induced refractive index change is no longer confined to the photosensitive region of the fiber 
(as is the case for FBGs written with ultra-violet lasers), but is determined solely by the fo- 
cusing geometry of the ultrashort laser. Both the size and position of the modifications depend 
heavily on the alignment and focusing of the writing beam [11,12]. Highly localized FBG ex- 
hibit a pronounced coupling to cladding modes. Our aim is to investigate and understand these 
cladding mode spectra in detail in terms of a full vectorial coupled mode theory. Our work 
shows that excitation of cladding-modes with strong selectivity and precision is possible, de- 
tailing how the coupling behavior is determined by the fiber geometry and to which extend it 
can be manipulated with the cross-section of the FBG. 
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2. Realization and spectra of highly localized FBG 

2.1. Fabrication and characterization 

The details of the grating inscription procedure closely follow that, which has been reported 
in [1 1, 13]. A femtosecond laser (Spectra-Physics, Hurricane) that delivered pulses shorter than 
110 fs at a wavelength of 800 nm was utilized for grating inscription. A 20 x oil immersion 
objective (NA = 0.8) was used for focusing the ultrashort pulses into the core of the fiber, as 
well as for imaging the fiber core before and after inscription of the gratings. All experiments 
were performed with a standard single mode fiber (Corning SMF-28e). The fiber was com- 
pletely stripped of its polymer cladding prior to grating inscription and subsequently threaded 
through a glass ferrule which was placed in close proximity to the objective lens. The lens, 
ferrule and fiber were immersed in an index matching oil in order to minimize refraction as 
the pulse propagates towards the focal spot. The ferrule was positioned such that the core of 
the fiber was located at the focus of the laser beam and the fiber was then pulled through the 
ferrule at a constant velocity while the grating was inscribed. By adjusting the ferrule's location 
with respect to the focus of the lens, the transverse position of the modifications of the FBG 
could be controlled to micrometer accuracy within the fiber core. The error in the longitudinal 
position of the modifications was much smaller, because the fiber was pulled by an air-bearing 
translation stage designed for constant velocity traversals, and the repetition rate of the laser 
was constant. All gratings were inscribed in SMF-28e fibers and are 20 mm long. Except where 
stated otherwise, the gratings have a period A = 1.062 um that reflects 1.5 um light in the 
second order. 




Fig. 1. Transmission differential interference contrast microscope images of the index mod- 
ifications viewed from (a) above (parallel to the direction of the writing laser beam) and 
(b) the side (perpendicular to the direction of the writing laser beam), (c) Schematic di- 
agram of the cross-section of the femtosecond induced modifications within the core of 
radius a, , establishing the coordinate system, where z points along the fiber and x and y are 
perpendicular to the fiber axis. 

Pulse energies between 200 and 275 nJ were used, which provided enough irradiance within 
the focal volume to cause a micro-explosion in the fiber, leaving behind a micro-void and hence 
forming a type II-IR Point-by-Point (PbP) FBG [14]. The micro-voids are surrounded by a 
compressed region of increased refractive index. In previous work, we have demonstrated that 
the net effective refractive index change is negative, although this is not significant in the present 
work [13]. Microscope images of the grating modifications taken from the top (a) and side (b) 
views are shown in Fig. 1. The micro-void regions, which are highlighted in Fig. 1, measure 
0.4 um in width by 1 .9 um in height. They are enclosed by a densified Type I shell of width 
approximately 1 um and height 8 um, where the index is slightly increased. The modifications 
are not necessarily centered, thus having an offsets of d x and d y . 



3 



2.2. Spectral response 

The measurement of the transmission and reflection properties was carried out with a swept 
wavelength system (SWS) (JDS Uniphase), that consisted of a narrow line-width external cav- 
ity laser diode. The probing wavelength of the laser diode can be continuously swept from 
A = 1520 to 1570 nm with a resolution of 3 pm. A typical transmission spectrum is shown in 
Fig. 2. In contrast to FBGs written with conventional techniques, the cladding mode structure 
is very pronounced: the reflectivity of core light into cladding modes is higher than 70 percent, 
even at wavelengths 20 nm shorter than the Bragg wavelength Xg = 1540 nm. For gratings at 
1555 nm, we see strong resonances 35 nm below the Bragg peak, limited only by the wave- 
length range of our measurement system. Comparable cladding mode spectra can be obtained 
with conventional techniques only with blazed FBGs with large tilt angles, e.g. 80° [6], whereas 
our gratings are untilted. 

As detailed in section 3, each resonance can be attributed to the coupling of the core mode to 
one or more definite cladding modes characterized by their radial and azimuthal mode numbers. 
It can be seen that the spectrum contains three sets of resonances with distinct well-defined 
envelopes. We will show that the envelopes correspond to modes with even (green dots) or odd 
(red dots) azimuthal mode numbers. 

Wavelength (nm) 



20 1530 1540 1550 1560 15 




1521 1522 1523 1524 1542 1547 1552 



Wavelength (nm) 

Fig. 2. Transmission spectrum of a PbP grating operating in second order at 1540 nm and 
showing strong cladding mode resonances. According to the notation laid out in section 3, 
the colored dots indicate the envelopes of resonances with HE / = 1 (red), EH / = 1 (blue) 
and HE 1 = 2 (green).The lower graphs show expanded scales to illustrate the HE/EH 
splitting of modes with high radial order m (lower left), and the existence of first-order 
cladding mode peaks right up to the second-order Bragg peak (lower right). 

Close inspection reveals that many of the resonances, especially at wavelengths lower than 
1530 nm (see lower left panel of Fig. 2), actually consist of doublets (marked by red and blue 
dots respectively). No such splitting is observed at the low order (long wavelength) resonances. 
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These doublets result from the hybrid nature of the cladding guided modes, which will be 
discussed in detail in section 3. 

Finally, the lower right panel of Fig. 2 provides a striking illustration of the strength of the 
cladding mode coupling in these gratings. It shows the response at wavelengths longer than 
the second-order Bragg peak at 1540 nm. A sequence of regularly-spaced pairs of resonances 
1-2 dB deep, each themselves with a recurring splitting pattern, is clearly visible. We interpret 
these peaks as the tail of the first-order cladding mode spectrum, and they therefore correspond 
to resonances many hundreds of modes below the fundamental. The actual first order Bragg 
peak would occur in the mid-IR around 3 pm and is beyond our observable range due to losses 
and the spectral range of our near-IR characterization setup. 

In the following sections, we perform calculations to understand these types of cladding 
mode spectra and explain their features in precise detail. 

3. Classification and properties of modes in three layer fibers 

For the purpose of calculating and classifying the cladding modes we consider a standard, 
cylindrical step-index fiber supporting a single core-guided mode over the wavelength range 
of interest. The elliptical modifications introduce a small anisotropy to the core and the modes 
are technically birefringent, but this effect is insignificant for the cladding modes. We therefore 
treat the fiber as having perfect circular symmetry, with core radius a, and cladding radius a 2 . 
The refractive indices are n u n 2 , respectively. 

Exploiting the small difference in n, and « 2 , cladding modes are most often calculated with a 
two-layer model consisting of a solid rod of index n 2 and radius a, surrounded by an air cladding 
with index 1 . They are further often treated in the approximation of weak-guidance or linearly 
polarized (LP) modes [15]. This approximation is usually sufficient for computing the position 
of the resonances in the spectrum [16, 17], since the propagation constants of the true HE and 
EH hybrid modes of the step index fiber differ only by a very small amount Aj3 [15]. 

It has previously been found, however, that in order to accurately calculate coupling coeffi- 
cients between core and cladding modes, both these approximations must be avoided [1]. The 
core-cladding boundary influences the electric field distribution of the lower order modes espe- 
cially. Hence it is necessary to solve the cladding modes for the full three-layer structure [18]. 
Moreover, because fiber modes of higher azimuthal order are not linearly polarized, we need 
a full vectorial approach [15] to accurately explain splittings in the cladding-mode spectra that 
do not appear in the LP picture. 

In the remainder of the paper, the core radius a\ = 4.15 pm and the cladding radius «2 = 
62.5 pm. The refractive indices of the core, cladding and surrounding air are n, = 1.4670,« 2 = 
1 .4618 and «3 = 1 .0 respectively. Modes travel in the z direction with the propagation constant 
j8 = 2Kh/X, where h is the effective index. A mode is guided within the core if h > « 2 , and in 
the cladding if n 2 > n> «3. To excite the system, we launch a single linearly polarized HE^i 
core mode. 

3.1. Field expressions for the modes 

In our analysis, we make extensive use of the three-layer fiber model in the form laid out by 
Erdogan [1]. He adapted the full analytical solution of the three-layer treatment of Tsao [18], 
which yields all hybrid modes with the azimuthal and radial integer indices / and m. However, 
in [1], Erdogan focused solely on cladding mode reflections of conventional fiber gratings, for 
which only / = 1 resonances occur. In our work, higher order modes with / > 1 are needed. We 
use the notation employed by Erdogan, but express the azimuthal dependence in trigonometric 
rather than exponential form. In cylindrical coordinates (r, (j), z), the electric E and magnetic 
H fields of the cladding modes inside the core (r < a,) can be expressed in terms of Bessel 
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with the transverse wavevector u, = (2n / X) yjn 2 — h 2 . The constant Zq = >J jUq /eo ~ 376. 7i2 
is the electromagnetic impedance in vacuum. Note that in these expressions «,j3,M, and P all 
depend on the mode indices / and m, and we sometimes indicate this explicitly by including the 
mode indices as subscripts. In particular, the mode parameter 



(2) 



characterizes the relative strength of the longitudinal field components, and is used to classify 
modes as HE or EH. The imaginary parameter £o is determined by solving the dispersion equa- 
tion for the three-layer fiber. The normalization constant Ef m is set such that all modes carry the 
same power. Complete expressions for the dispersion relation and the cladding mode fields in 
the regions r > a, are provided in appendices A and B respectively. 

Note that Eqs. (1) represent two orthogonal sets of solutions, distinguished by the rotation 
angle <p, which we take as <p = or <p = — n/2. Thus, all hybrid mode solutions of the form 
in Eqs. (1) appear as degenerate pairs of fields with orthogonal polarization states. Since in our 
experiments the input mode HEi i is y-polarized, it is represented purely by the q> = state. In 
general, a grating may produce coupling to cladding modes of both q> = and <p = — n/2. As 
discussed below however, careful alignment of the polarization as well as the defect orientation 
and location eliminates coupling to the perpendicular (p = —n/2 modes. 

3.2. Classification of the modes 

It is useful to recall some standard properties of the core-guided hybrid modes of the two-layer 
model. In that system, the modes are classified as HE or EH according to the value of P[ m . The 
HE modes have — 1 < P/„, < and the EH modes have P;„, > 1, with P/ m m ±1 for the strongly 
bound modes, (and P\ m — ±1 in the LP limit). For azimuthal number /, the transverse field 
components of the HE modes have predominantly a 7/_i(«,r) radial dependence, whereas the 
EH modes are dominated by the («,r) Bessel function. Consequently, a mode HE/ „, has a 
time-averaged Poynting flux (S z ) with peaks closer to the origin than the corresponding EH/ ,„ 
mode. For fixed /, the modes occur in a strict alternating sequence of HE and EH, (excepting 
the fundamental HEi.i . 

The mode parameter P/,„ also turns out to be very convenient for the classification of cladding 
modes in the three-layer fiber. However, the HE/EH criterion has to be modified. For further 
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discussion, we define the projected normalized wavevector 



which corresponds to the normalized frequency a cladding mode would have in the core region 
if the refraction on the core-cladding boundary is ignored. 

In Fig. 3, we plot the hybrid mode parameter P/„, against the normalized wavevector for the 
lowest six azimuthal orders. The plots show that P\ m has a complex dependence and is not 
restricted to the ranges — 1 < P < and P > 1. Adjusting the two-layer fiber definitions, we 
designate the modes with \P/ m | < 1 as HE (colored blue) and those with \P[ m | > 1 as EH (colored 
green). Note that with this definition, we retain consistency with the classification of core modes 
in the two-layer picture, and the modes again appear in a strict sequence alternating between 
HE and EH. Now in Fig. 3 we find that for each value of /, the diagrams can be divided into 
regions according to the sign of P[ m for either set of modes. In fact, the transitions, indicated 
by vertical dashed lines, correspond to mode indices to at which an additional intensity ring 
enters the core. We refer to these discontinuities as "virtual cut-offs" t//,,/, where the virtual 
radial number to' is to be associated with the number of maxima within the core. The cutoff 
value U[ m i is the value of Ui m at which the m'th maximum appears in the core. The intensity 
distributions near and in the fiber core S z = \ (E r H^ — E^H*) for a number of modes are plotted 
as insets. It can indeed be seen that a further intensity ring appears within the core each time 
there is a discontinuity in the P\ m plot. 

The character of the modes changes qualitatively at the first few virtual cutoffs. Cladding 
modes with Ui,„ < 0/ j/ have essentially no intensity within the core. (As with the classic two 
layer model, U\y = 0. Therefore, / = 1 modes always guide some light within the core region.) 
When the normalized wave-vector satisfies U[ \i < U < U[ 2', HE modes {P[ m « — 1) have one 
maximum within the core. In contrast, EH modes (P/ m « 1) have a minimum at the center 
and have significant field contributions only near the core-cladding boundary. Thus the spatial 
distribution of energy is distinct for HE versus EH. The polarization of each mode in this regime 
is quite complex and far from uniform across the fiber. 

The propagation of modes with U > 0/ 2 is more determined by total internal reflection 
at the cladding/surrounding interface than the core/cladding interface. In this region we have 
either \Pi,„\ ~ (HE), or \Pt,„\ ^> 1 (EH), which is indicative of modes for which the energy 
in the longitudinal components resides almost purely in H z (HE) or purely in E z (EH). As one 
or other longitudinal component is close to zero, consequently the transverse field of the HE 
modes eventually become close to purely azimuthally polarized (quasi-TE) and that of the EH 
modes becomes strongly radially polarized (quasi-TM) [19]. 

4. Calculation of the spectra 

We now turn to the detailed explanation of the form of the cladding mode spectra, calculat- 
ing the resonant wavelengths, the coupling coefficients and finally the complete transmission 
spectra. 

4.1. Resonance wavelengths 

One of the challenges of computing FBG spectra of strong cladding mode coupling is that 
it usually occurs over large bandwidth ranging from 20 to several hundred nanometers. All 
resonances that result from the reflection of the incoming mode HEi 1 into cladding modes 
(7, to) by the grating of period A are defined by the phase matching condition 

/Jn±/J /m -27r/A = 0, (4) 
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Fig. 3. Hybrid mode parameter P\ m for I = 1 to / = 6 as a function of the projected normal- 
ized wavevector (//,„. Blue and green circles denote modes labeled HE and EH respectively. 
The insets display the intensity of the cladding modes in the region r < 3a i with the black 
circle denoting the core boundary (EH modes are shown above the x-axis, HE modes below 
it). Vertical lines indicate the virtual cutoffs U[ m i and are labelled on the upper edge of each 
graph. The change in character of the modes at each virtual cutoffs is clear. 



Taking all the propagation constants as positive, the upper sign holds for FBGs, and the lower 
sign for long period gratings. 

Since Eq. (4) involves the wavelength implicitly through the effective indices of the modes, 
computing the resonance wavelengths requires solving the dispersion relation for hundreds of 
modes at finely discretized wavelengths. Within a bandwidth of 30 nm, the material dispersion 
of the fiber can usually be neglected as it is of the order of 10~ 4 and affects both core and 
cladding refractive indices in a similar fashion. Therefore the index contrast n, — «, varies only 
weakly. However, the waveguide dispersion is significant and must be accounted for. Solving 
Eq. (4) exactly is relatively straightforward, but considerable care is required in solving the 
dispersion relation (11) to find all the required effective indices h\ m at all wavelengths, since 
some of the EH modes appear as resonances as narrow as one part in a million. 

If the intention is to identify particular modes in a measured spectrum rather than precisely 
predict their location in advance, a more convenient approximate approach that provides addi- 
tional insight is to solve the dispersion relation at the Bragg wavelength Ao and extrapolate the 
waveguide dispersion of the propagation constants using the relation [20] 

_ ... ( 2 (2tz\ 2 V /2 (/+2/n+l/2)V 

ftm(A) ~ y n i - A imj with Ai m = ± 4 ^ , (5) 

where A/ m is a mode-dependent dispersion-free constant. This relation applies for the two-layer 
model of a fiber in air («3 = 1) in the LP approximation and the regime of high V number 
(for the cladding modes V « 270). These LP resonances represent the HE/EH doublets for 
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resonances of high radial number m described earlier. Although this approximation neglects 
the vectorial mode splitting and is rather inaccurate for low m, it is helpful for understanding 
the degeneracies with respect to the azimuthal number /: because of the / + 2m term, two sets 
of propagation constants exist, one for odd I and one for even / halfway between. 

In order to precisely predict the resonance wavelengths using Eq. (5), for each mode the 
dispersion free constant 

, _ 2 (2n\ 2 

Mm = («2 -«/,„) It:)' (6) 

was computed by solving the exact three-layer model for «/, n (Ao) at just the Bragg wavelength 
Aq. Substituting A/ m for A/ m in (5) and combining that equation with (4), we obtain a quadratic 
equation for the wavelength of the (l,m) resonance 



nil _i_ 
A 31 








1/2 









(7) 



The positive root gives the wavelength for coupling into counter-propagating modes. The neg- 
ative root yields the resonances for coupling into co-propagating modes and is therefore ne- 
glected in the following treatment. We have checked that this extrapolation procedure provides 
results in very close agreement with the exact approach. Figure 4 depicts a similar transmis- 
sion spectrum to that in Fig. 2, together with the three-layer multimode resonances, computed 
according to the steps just described. The modes were calculated at Ao = 1555 nm. As ex- 
pected from the LP approximation in Eq. (5), the numerical solution also yields two sets of 
nearly-degenerate resonance doublets, one for odd / (red lines) and one for even / (blue lines). 
For given /, each doublet corresponds to an HE/EH pair. The alignment of the measured and 
calculated resonances is good. For low wavelengths (Fig. 4(b)), there is a small wavelength 
shift compared to the measured spectra which arises because of the neglected material disper- 
sion, however the splitting of the HE/EH doublets is in very close agreement. In this range, 
the small difference in the propagation constant Phe — Peh of the HE/EH doublets is con- 
stant, as expected from the standard 2-layer hybrid mode solution applied to the cladding-air 
boundary [15]. This is not the case for cladding resonances closer to the Bragg peak (Fig. 4(c)): 
Phe — Peh becomes very small in the vicinity of a virtual cut-off. It can also be seen that in 
this regime, the doublets of a different azimuthal order are shifted such that the EH resonance 
of a lower I is in most but not all cases degenerate to the HE resonance for I + 2. The use 
of LP modes LPf m = HE^ + i +EH^_j „, still leads to surprisingly good results for coupling 
strength of lower m cladding modes [17], since the contribution of higher / modes cancels out 
in transversally homogeneous FBGs as shown in the following sections. 

4.2. Computation of coupling constants and spectral response 

The interaction of the modes with the FBG can be computed by solving the coupled mode 
equations as established by Kogelnik [21]. The unperturbed mode fields of the fiber and their 
propagation constants have been calculated beforehand as described in the preceding sections. 
Here, we use the coupled mode equations adapted for FBGs by Erdogan [1]. This set of differ- 
ential equations is determined by the coupling constants, which are obtained by evaluating an 
overlap integral of the interacting modes with the perturbation. In the case of FBGs, the trans- 
verse coupling constant that governs the coupling of the incoming fundamental mode HEi i to 
any other reflected mode (l,m) is computed with 

4 CO f 271 I"" T T* 

K im = ~ * T • / d(j> drrA £ (rJ)Ej r Ej m , (8) 
nv 4 Jo Jo 



9 




Fig. 4. Measured cladding mode spectra in transmission (bottom) and computed HE and EH 
resonances (top) for odd (red) and even (blue) azimuthal order I. (a) Complete measured 
wavelength range; (b) low-and (c) high-wavelength range in more detail. 



where only the transverse electric fields E /m are evaluated. (The coupling associated with the 
longitudinal field components is orders of magnitude smaller and can be neglected.) The per- 
turbation of the dielectric constant Ae(r, <p) is induced here by refractive index modifications of 
the grating. The grating is modeled as periodic but non-sinusoidal in the propagation direction. 
For simplicity, a rectangular profile with a 50% duty cycle is assumed, and the factor 4/ (nv) 
accounts for the Fourier component at reflection order v. The investigated FBG reflects in sec- 
ond order (v = 2). To alleviate the need to compute a system of several hundred coupled mode 
equations (one for each cladding mode), the synchronous approximation is used in which all 
modes not near resonance are neglected. Solving the remaining equations yields wavelength 
dependent amplitudes, from which reflection and transmission spectra can be accurately com- 
puted [1,4]. 

4.3. Fundamental coupling properties of the modes 

Before characterizing the modification term Ae in Eq. (8), it is useful to elaborate on the kernel 
K = Ae(r, 0)Ejj ■ Kj m of the coupling constant integral, since this can be evaluated before inte- 
grating with the help of the results of section 3. Here, the incoming mode is always the HEu 
core mode with <p = 0. This mode couples with any other other hybrid mode, characterized by 
P[ m and U[ m = {2%/X)yJn} — h 2 with 

K{<p) = Ae(r 7 (j))EnEi m L -^j^-Jo(u u r){(l -P, m )//_i (« /m r)cos[(/ - 1)0 + <p] 

-(1 +Pi m )Ji +i (u, m r)cos[(l +l)^ + (p]} (9) 

as the kernel. Note here, that the core mode can couple to modes of both sets of polarization 
orientation, determined by <p. We refer to the parallel kernel = K((p = 0) for coupling to 
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modes with the same azimuthal dependence as the launched core mode, and the perpendicular 
kernel K^- = K(q> = —iz/2) for coupling to the orthogonal set of modes. 

If the FBG extends homogeneously over the whole cross section of the fiber core and its 
surrounding, Eq. (8) would be equivalent to the orthogonality relation of the modes, evaluating 
to zero for any mode (l,m) other than the core mode HEi j. In general, coupling between 
modes of different indices occurs if these conditions are violated. In conventional FBGs for 
example, the grating is only formed in the photosensitive core and not the cladding. This case 
corresponds to a trivial integration over <p and any terms with a sine or cosine cancel out. Thus, 
only the first term of Eq. (9) remains for I = 1 and arbitrary m, while the perpendicular kernel 
K L always sums to zero. This case is described in detail in [1]. In particular, coupling to higher 
order modes with I > 1 is only possible if the refractive index of the FBG is inhomogeneous 
or asymmetric. The strength of such coupling is influenced by the cross section of the FBG 
through its non-uniformity and localization. 

The classification of the fiber modes as discussed in section 3.2 is reflected in the kernel, 
too. For the cladding modes below the U/ 2' cutoff, the strength of the coupling coefficient has a 
pronounced alternating character as a function of m due to the hybrid mode parameter P/,„ « ±1 : 
The (1 —Pim) term in Eq. (9) is responsible for the strong selective character. While for HE 
modes the J[_\(ui m r) term dominates, it is the // + i(w/ m r) term for EH. In contrast, cladding 
modes beyond the XJiji cut-off have quasi-TM/TE core fields. Because of this, P/„, ^ ±1 and 
HE and EH modes differ less in their coupling behavior. 

4.4. Coupling coefficients and spectra 

We now proceed to calculate the coupling coefficients and complete transmission spectra for 
highly-localized FBGs. In our experiments, the FBG consists of femtosecond pulse induced 
micro voids (see Fig. 1 and Ref. [13],) that are approximately prolate spheroids (with semi- 
major axis 0.95 um and semi-minor axis 0.2 |jm,). In the following, the densified shell around 
the microvoid is neglected, since its influence on the overall coupling is small compared to the 
void itself [13]. To reduce the geometry to the form of Eq. (8), we treat the grating as a periodic 
square modulation with an elliptical cross-section and an index contrast obtained by averaging 
the local index of the true structure over a whole period. Assuming this, the perturbation of the 
dielectric constant is 

Ae(x,y) = £o(«L/ - n\)0{x,y) (10) 

where 0(x,y) = 1 inside the ellipse and elsewhere, and n moc ( = 1.325 [13]. 

In our fabricated gratings we can position the voids arbitrarily in the core (although there is 
some variation in the void displacement along the core due to vibration in the ferrule during 
the writing process). The strength of the cladding mode coupling strongly depends on small 
displacements of the micro voids from the central axis of the fiber (results presented below 
in Fig. 6). For inhomogeneous FBGs, the launch polarization has a significant influence on the 
overall coupling strength because of coupling to both sets of modes (Eq. (9)). These polarization 
dependent effects will be discussed in more detail in a forthcoming publication. In the following 
calculations, we emphasize more the dependence of the coupling to the radial position of the 
highly localized FBG. We therefore choose the special case for which the launch polarization is 
aligned with the long axis of the modifications and the modifications are only displaced along 
the x-axis. The perpendicular kernels K therefore cancel out and can be neglected. 

The coupling integral was evaluated for displacements d x in the jc-direction, ranging from 
to the core cladding boundary at 4.15 \im. The resulting coupling constants for the parallel 
kernel K" are plotted in Fig. 5(a) as a function of radial mode order m and displacement d x for 
azimuthal orders / =1-6. The horizontal white dotted lines show the locations of the virtual 
cutoffs for each value of I. Since resonances above JJq y lie beyond the wavelength range of our 



11 



(a) Azimuthal index (I) 
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Fig. 5. (a) Core-cladding mode coupling coefficients for different displacements d x of the 
microvoid in x-direction from the center of the fiber core. Coupling constants for / = 1 
to / = 6 are plotted, m is ascending from top to bottom. Dotted lines indicate the virtual 
cutoffs, (b) Computed transmission spectra for different displacement of the modification 
in the .^-direction. 



measured spectra, we expect negligible coupling to modes of / > 6, so that the range / =1-6 
is appropriate. All of the concepts of the preceding section are present in these plots: there is 
almost no coupling below the C/j y cutoff, pure HE coupling for modes below U[ y (as indicated 
by the bright red bands for odd m (HE) and dark bands for the even m (EH) modes); and reduced 
contrast between the HE and EH selective coupling above £7/2', when the nature of the modes 
changes to quasi-TE and quasi-TM. In agreement with [1] and the earlier symmetry argument, 
there is no coupling into modes of / ^ 1 if the modification is centered (d x = 0). Figure 5(b) 
displays the corresponding spectra calculated from the coupling constants of 5(a) as outlined in 
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section 4.2. The even / resonances appear in the spectra in between the odd / resonances (see 
Eq. 5). Most striking is the sensitivity of coupling to / = 2 modes: even a small displacement 
d x ~ 0.4 |jm of the modification suffices for measurable coupling strength. The transmission 
spectra do not allow us to distinguish the strength of coupling to modes with I of the same 
parity, but the maps of coupling coefficient in (Fig 5 a)) indicate that measurable coupling to 
/ = 3 and / = 4 modes can only be observed if the modification is at least 1 iim from the center. 



(a) 




1536 1538 1540 1542 1544 1546 1548 1550 1552 1554 1556 

Wavelength [nm] 



Fig. 6. Measured (a) and computed (b) transmission spectra for the micro void at d x = 0.8 
d y = 0.7. The insets show the EH/HE doublets. 



5. Discussion 

In Figure 6 we compare measured and calculated transmission spectra for one FBG with the 
light launched perpendicular polarized to the fiber axis. The positions of the FBGs within the 
fiber core were estimated from microscope images. Due to fiber vibration in the ferrule, the 
micro voids are not uniformly distributed along the FBG, and the displacement values used in 
calculations are the maximum measured values. In the first case (Fig. 6(a) and (b), the spec- 
trum was simulated for d x = 0.8 |jm and d y = 0.7 iim. Both in experiment and simulation, the 
polarization direction was along the long axis of the micro void, parallel to the y- axis. Because 
of the displacement in y-direction, the perpendicular coupling constants (<p = — n/2 in Eq. 9) 
had to be considered as well. The computed spectra closely agree with the experimentally ob- 
tained spectra. The slight difference in the envelope can be attributed to the restriction of the 
model to account only for uniform displacement. (Note that we have also performed coupling 
calculations using a core mode that was calculated using the finite element method to include 
the mean index perturbation associated with the grating modifications. This makes a negligible 
change to the transmission spectra.) 

Due to the displacement from the center and the therefore asymmetric cross-section of the 
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FBG, coupling to higher order modes occurs. This is clearly visible in the second comb of res- 
onances in between the / = 1 resonances, which corresponds to the even / modes [see Eq. (5)]. 
Additionally, with the aid of imaging the reflected cladding modes, we were able to determine 
that interference between degenerate modes of different azimuthal order took place and there- 
fore that higher order azimuthal modes were significantly excited. This will be the subject of a 
future publication. 

Consistent with the results in Fig. 4, the low order resonances in Fig. 6(a) appear as single 
peaks but doublets appear at higher m (shorter wavelengths). The single peaks correspond to 
pure HE modes. The EH modes are formally present in these peaks, but as seen in Fig. 5(a), 
the coupling is too weak for them to be observed. While the single peaks show almost no po- 
larization dependence, the strength of the doublet peaks heavily depends on how the axis of 
polarization is orientated with respect to FBG. The experimental and theoretical results con- 
cerning this aspect will also be presented in a future publication. 

6. Conclusion 

In this paper we investigated the core-cladding coupling behavior of FBGs that have trans- 
versely inhomogeneous refractive index profiles. For that purpose, we wrote intra core FBGs 
in step index fibers based on the point-by-point approach. The micro voids formed within the 
index modifications that the grating comprises, are much smaller than the core and therefore 
allow for an accurate excitation of particular modes. Here, even I resonances served as a strong 
test for the centeredness of a FBG. Our findings were validated by a rigorous calculation of 
the coupling constants. Accurate spectra of transversely inhomogeneous FBGs could be re- 
produced, which predicted many of the transmission characteristics of such FBGs, that have 
only been described qualitatively before. The necessary theoretical foundation for a detailed 
investigation of polarization dependent features could be established. 

We have demonstrated efficient mode converters to odd and even I modes, realized with 
off-centered, highly localized FBGs. Such intra-core gratings allow for specifically targeting 
coupling to cladding modes of higher azimuthal index /. We believe these findings will enable 
many new fiber sensor designs, that rely on asymmetrically structured fiber cores. The efficient 
core-cladding mode conversion presented here can also be exploited for in-fiber elements that 
allow for the tailored excitation of higher-order mode fibers [22]. Our extensive treatment can 
be easily transferred to fibers that carry multiple modes in the core for design and calculation 
of core-mode converting gratings. 
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A. Dispersion relation 

According to [1,18], the dispersion for the three layer fiber with arbitrary integer azimuthal 
number / is 
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Here, we have defined 

uf m = {2n/X) 2 (n 2 ~h 2 ), u\ = {2jc/X) 2 {n 2 -n 2 ), w 2 = (2n/X) 2 (h 2 ~ n 2 ), (12) 
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The dispersion relation includes Bessel-functions of the first and second kind J„ and N„ as well 
as modified Bessel-functions K n of the second kind. In addition, the products 

pi{r) =Ji{u 2 r)Ni{u 1 a l )—Ji{u^ti)Ni{u 2 r), (14a) 

qi(r) = J/^^NKu.a,) ~ j',(u 2 ai)Ni(u 2 r), (14b) 

ri{r) =j' l (u 2 r)Ni{u 2 a l )-Ji{u 2 a l )N',(u 2 r), (14c) 

si{r)=J' l {u 2 r)N' l {u 2 a 1 )-J' l (u 2 a 1 )N[{u 2 r), (14d) 

are used, where the prime stands for differentiation with respect to the total argument. 
B. Fields 

Here we provide the expressions for the fields in the cladding and air regions, Inside the 
cladding (a, < r < a 2 ): 
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In the air region, (r > a 2 ) the fields are 

Ef=-E lm C lm — ^^ ^( Ws r)<dn(/» + y) e ^-») 



,2 



11 



iEj m Ci m 



Ki(w^a 2 )w 3 
Ki(w 3 a 2 )w 3 



l3-K,(w s r) + ^Kl(w 3 r) 
w 3 r Inf 

w 3 rnf 



cos(Z0 + ^)e'' ( ^-' B ' ) 



(18a) 
(18b) 
(18c) 



H 



h: 

with 



cl 



iE\ m C\, n 



<JFt,W 3 



ZqKi(w 3 ci 2 )w 3 Pi 



- cos 



11 



-El m Q 



l 



Zq K[(w 3 a 2 )w 3 
1 u\ 
ZQKi(w,a 2 )w 3 



IGt. , , Ft.o ., 
rw 3 / 

-G 3 ^'(w 3 r) + — X,(w 3 r) 



(19a) 

20s{l$ + (p)e i[fiz - m) (19b) 
sin(/(/) + ^) e , ' ( ^-' a ' ) (19c) 



F 3 



-F 2 pi(a 2 ) + —qi(a 2 



GlPl{a 2 ) ?—qi{a 2 



(20) 



16 



